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Abstract 

In this paper, we present a model of the formation and chemical enrichment of 
I elliptical galaxies that differs from the conventional picture in two fundamental ways: 

'sj" ■ 1. Ellipticals do not form in a single monolithic collapse and burst of star formation 

I at high redshift. Instead, most of their stars form at modest rates in disk galaxies, 

' which then merge to form the ellipticals. 

I 2. Galaxies do not undergo "closed-box" chemical evolution. Instead, metals can be 

' transferred between the stars, cold gas and the hot gas halos of the galaxies. It is 

Q_|. assumed that metals are ejected out of disk galaxies during supernova explosions 

Q I and these metals enter the hot gas component. The fact that metals are more 

^ ' easily ejected from small galaxies leads to the establishment of a mass-metallicity 

c/3 . relation for the disk systems. Big ellipticals are more metal rich because they are 

^ ' formed from the mergers of bigger disks. 

- We use semi- analytic techniques to follow the formation, evolution, and chemical enrich- 

r> I ment of cluster elliptical galaxies in a merging hierarchy of dark matter halos. The in- 

' elusion of the new metallicity-dependent spectral synthesis models of Bruzual & Chariot 

enable us to compute the colours, line indices and mass-to-light ratios of these galaxies. 
We find that with physically realistic parameters and with the assumption that feedback 
is efficient, even in massive galaxies, we are able to reproduce the slope and scatter of the 
colour-magnitude and the Mg2 — a relations. We are not able to match the increase in 
M/L for bright ellipticals if we assume their mass to measure purely the total quantity 
of stars in the galaxy. We also study the evolution of these relations to high redshift. 
We show that the luminosity- metallicity relation does not change with redshift, but the 
mean stellar age of the galaxies scales with the age of the Universe. This is why the 
evolution of cluster ellipticals appears to be well described by simple passive evolution. 
Finally, we study the enrichment history of the intracluster gas. Our models predict that 
more than 80% of the metals were ejected by galaxies with circular velocities less than 
250 km s~^ at redshifts greater than 1. The metallicity of the ICM is thus predicted to 
evolve very little out to z > 1. 
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1 Introduction 



The relationship between the colours, line strengths and luminosities of elliptical galaxies is 
potentially a source of information about physical processes important during their formation. 
The conventional interpretation of the colour-magnitude relation of ellipticals is that it is 
primarily a metallicity effect (Faber 1977; Dressier 1984; Vadcr 1986). It is often assumed 
that elliptical galaxies form monolithically in a single giant burst of star formation at high 
redshift. Supernova explosions heat the gas, and when its thermal energy exceeds its binding 
energy, a galactic wind is produced, which is assumed to interrupt further star formation 
and enrichment. Since the binding energy per unit mass of gas is higher in more massive 
galaxies, these are able to retain their gas for much longer times and so to reach higher 
metallicities than less massive galaxies (Larson 1974; Matteuci & Tornambe 1987; Arimoto & 
Yoshii 1987; Bressan et al 1996). This interpretation is not unique, however, because colours 
and line strengths are degenerate in age and metallicity. As a result, the colour-magnitude 
relation of local ellipticals can also be interpreted as an increase in mean stellar age with 
luminosity (Worthey 1996). In recent work, Kodama & Arimoto (1997) have pointed out that 
the evolution of the colour-magnitude relation to high redshift does not support the hypothesis 
that age effects alone arc responsible for the redder colours of the brighter ellipticals. Following 
this. Chariot, Ferreras & Silk (1997, in preparation) have explored models where both age 
and metallicity are important, and have shown that several of these can provide equally good 
fits to the data. 

An alternative approach to interpreting the spectra of ellipticals is to begin with a more 
realistic dynamical picture for the formation of these objects, and then to study the properties 
of their stellar populations by applying spectral synthesis techniques. Recently, there has been 
considerable progress in modelling the observed properties of ellipticals in a scenario where 
they are formed by the merging of disk galaxies in a universe where structure is built through 
hierarchical clustering (Kauffmann 1996b, Baugh, Cole & Frenk 1996). This work has, for the 
first time, placed the theory of the formation and evolution of elliptical galaxies in a proper 
cosmological context. In these models, star formation takes place mainly in centrifugally- 
supported disks, formed when gas cools and condenses at the centres of virialized halos of dark 
matter. A self-regulated balance between cooling and feedback from supernova explosions 
maintains star formation at a roughly constant rate in these disks. If two disk galaxies 
merge once their surrounding halos have coalesced, an elliptical galaxy is formed and all cold 
gas present in the merger remnant is transformed instantaneously into stars in a "burst". 
Kauffmann (1996b) showed that although ellipticals are predicted to have been assembled by 
relatively recent mergers {z < 1) in a standard cold dark matter (CDM) cosmology, most of 
their stars form at considerably higher redshift in their disk-galaxy progenitors. The K-band 
luminosity-weighted stellar age of cluster ellipticals was shown to vary between 8 and 12 Cyears 
and the scatter in the predicted V — K and U — V colours was of order 0.04, consistent with the 
small observed scatter of Virgo and Coma cluster elliptical colours (Bower, Lucey & EUis 1992). 
The model presented in that paper did not include any scheme for the chemical enrichment of 
galaxies and it was unable to reproduce the slope of the observed colour-magnitude relation. 
An important conclusion of this work was that age effects in hierarchical models could not 
alone account for the redder colours of the more massive ellipticals. 

We have now included a simple prescription for chemical evolution in our semi-analytic 
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galaxy formation scheme in order to analyze the effect of variations in metallicity on the 
observed properties of ellipticals. Supernovac take place in stars forming disk galaxies and 
the explosions cause cold gas and metals to be ejected into the halo, where the metals are 
mixed with the hot halo gas. The efficiency with which this process occurs depends on the 
potential well depth, so that less massive galaxies eject more material and hence produce 
fewer stars and metals per unit mass of available cold gas. The metallicity of disk galaxies 
thus increases as a function of their mass. We then show that in a hierarchical merger picture, 
bright ellipticals form from mergers of substantially more massive disks than faint ellipticals. 
This lays the theoretical groundwork for understanding the origin of a mass-metallicity relation 
for elliptical galaxies. 

We then compare the predictions of the model with observational data. To do this, we 
make use of the new stellar population synthesis models of Bruzual & Chariot (1997, in 
preparation); these enable us to compute the colours and spectral line indices of evolving 
stellar populations of different metallicities. We show that with physically reasonable choices 
for the free parameters in our model, we are able to reproduce both the slope and the scatter 
of the U — V and V — K colour-magnitude relations of cluster ellipticals. In addition, we obtain 
a good fit to the slope of observed relation between the Mg2 index and the stellar velocity 
dispersion a, where it is assumed that a simply scales in proportion to the circular velocity 
of the halo in which the elliptical galaxy formed. We also examine the evolution of these 
relations to high redshift and show that out to a redshift of 1 in a standard CDM cosmology, 
the slopes change very little and the colours and index values at fixed luminosity shift to lower 
values as expected for a passively evolving population of stars. At redshifts greater than one, 
there is a slight flattening in the slopes of the relations, caused by the slightly younger mean 
ages of the bright ellipticals relative to the faint ones. 

Although the colour-magnitude relations are reproduced very well by our model, we can- 
not account for the observed increase in the mass-to-hght ratios of bright elliptical galaxies if 
we assume the mass, estimated using the central line-of-sight velocity dispersion a, is purely 
a measure of the total quantity of stars in the galaxy. We hypothesize that since massive 
ellipticals are built up through a greater number of mergers than their less massive counter- 
parts, more dark matter may end up being mixed with their stars. Finally, since our chemical 
evolution model tracks the recycling of heavy elements between stars, cold gas, and the hot 
intergalactic medium, we are able to study the enrichment history of the hot gas in clusters. 
With a Sn la + Sn II yield of twice solar, we are able to reproduce the observed ratio of iron 
mass to total blue luminosity in nearby clusters (Renzini et al 1993). We show that 80% of the 
metals are produced at redshifts greater than 1 in galaxies with circular velocities less than 
250 km s~^. The metaUicity of rich clusters is predicted to evolve very little out to redshifts 
of 1 and beyond. 

2 Modelling the Formation and Chemical Enrichment 
of Elliptical Galaxies 
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2.1 Semi-analytic models of elliptical galaxy formation in a hierar- 
chical Universe 

The semi-analytic model we employ is the same as that used by Kauffmann (1996b) to compute 
the stellar ages of elliptical galaxies and bulges formed by mergers in a cold dark matter 
universe. The reader is referred to that paper for more background. The new material in the 
present paper is a prescription for chemical enrichment and the transfer of metals between 
stars, cold galactic gas and hot halo gas. We describe this in detail in the next section. 
More details about semi-analytic techniques may be found in Kauffmann & White (1993), 
Kauffmann, White & Guiderdoni (1993, KWG) , Cole et al (1994) and Baugh, Cole & Prenk 
(1996). 

To summarize: 

An algorithm based on an extension of the Press-Schechter theory due to Bower (1991) and 
to Bond et al (1991) is used to generate Monte Carlo reahzations of the merging paths of dark 
matter halos from high redshift until the present. This algorithm allows all the progenitors 
of a present-day object, such as a cluster, to be traced back to arbitrarily early times. Dark 
matter halos are modelled as truncated isothermal spheres and it is assumed that as the 
halo forms, the gas relaxes to a distribution that exactly parallels that of the dark matter. 
Gas then cools and condenses onto a central galaxy at the core of each halo. Star formation 
and feedback processes take place as described in KWG (see also below). In practice, star 
formation in central galaxies takes place at a roughly constant rate of a few solar masses per 
year, in agreement with the rates derived by Kennicutt (1983) for normal spiral galaxies. 

As time proceeds, a halo will merge with a number of others, forming a new halo of larger 
mass. All gas which has not already cooled is assumed to be shock heated to the virial 
temperature of this new halo. This hot gas then cools onto the central galaxy of the new 
halo, which is identified with the central galaxy of its largest progenitor. The central galaxies 
of the other progenitors become satellite galaxies, which are able to merge with the central 
galaxy on a dynamical friction timescale. If a merger takes place between two galaxies of 
roughly comparable mass, the merger remnant is labelled as an "elliptical" and all cold gas 
is transformed instantaneously into stars in a "starburst". Note that the infall of new gas 
onto satellite galaxies is not allowed, and star formation will continue in such objects only 
until their existing cold gas reservoirs are exhausted. Thus the epoch at which a galaxy is 
accreted by a larger halo delineates the transition between active star formation in the galaxy 
and passive evolution of its stellar population. The stellar populations of elliptical merger 
remnants in clusters hence redden as their stellar populations age. Central galaxy merger 
remnants in the "field" are able to accrete new gas in the form of a disk to form a "spiral" 
galaxy consisting of both a spheroidal bulge and a disk component. Note that as in KWG, 
we suppress the formation of visible stars in cooling flows within halos with Vc > 500 km s~^ 
in order to avoid producing central cluster galaxies that are too bright and too blue. 

The cosmological initial conditions used in this paper are a, as — 0.67 CDM universe 
with Q — 1 and Hq — 50 km s~^ Mpc~^. Although this model does not match the COBE 
measurements of the amplitude of the microwave background fluctuations, it has become the 
de facto standard for galaxy formation studies, since it provides a reasonable fit to observed 
small-scale galaxy clustering. 
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2.2 The chemical enrichment of galaxies in a hierarchical Universe 

In simple closed-box models of chemical evolution, all metals produced by stars are retained 
within the galaxy. After all the gas has been converted into stars, the metallicity of the 
galaxy is then simply given by the mass of metals produced per solar mass of long-lived stars 
that are formed. In our model, metals are exchanged between three components: stars, cold 
galactic gas, and hot halo gas. The routes through which this exchange may take place are 
the following: 

1. Stars — > hot and cold gas As decribed in KWG, the star formation rate in disk galaxies 
is given by 

M stars = <^Mcold/tdyn, (1) 

where Mcoid is the mass of cold gas in the galaxy, tdyn = Rvir/^^Vc is a measure of the 
dynamical time of the galaxy (this assumes that the galaxy collapses by a factor 10 from 
the virial radius of the halo), and a is a free parameter. We assume that stars form with 
a standard Scalo (1986) IMF truncated at 0.1 and 100 Mq and use the instantaneous 
recycling approximation (e.g. Tinsley 1980) to compute the amount of gas and metals 
returned to the interstellar medium and halo. A yield y of heavy elements is produced 
per solar mass of gas converted into stars. A fraction / of this yield is uniformly mixed 
with the hot gas in the halo, and the rest is uniformly mixed with the cold gas in the 
galaxy. The parameter / is introduced to allow for the possibility that heavy elements 
may be ejected directly into the halo in galactic fountains or chimneys resulting from 
the simultaneous explosions of tens or hundreds of supernovae (see Spitzer 1990 for a 
review). The gas fraction returned by evolved stars is estimated from the population 
synthesis models to be 25 % and is roughly independent of metallicity. This material is 
added to the cold gas component of the galaxy. 

2. Cold gas — > hot gas In our models, it is postulated that supernova explosions can 
release enough energy to drive cold gas into the halo. For a standard Scalo IMF, the 
number of supernovae expected per solar mass of stars formed is 7]sN IQ-^M-^. 
The kinetic energy of the ejecta from each supernova is 10^^ ergs. If a fraction e of 
this energy is used to reheat cold gas to the virial temperature of the halo, the amount 
of cold gas lost to the intergalactic medium in time At may be estimated from simple 
energy balance arguments as 

AM„,., = 4%^At (2) 

5 

Here e is a free parameter controlling the efficiency of the feedback process. When cold 
gas is returned to the halo, it is assumed that its metals are uniformly mixed with the 
halo gas. 

3. Hot gas — > cold gas Metals present in hot halo gas may incorporated into a galaxy 
when that gas is able to cool and condense at the centres of halos. The coohng rates 
in halos are calculated as described in KWG and are based on the model by White Sz 
Frenk (1991). The cooling functions vary according to the metallicity of the hot gas and 
are interpolated from Figure 9.9 of Binney & Tremaine (1987). 
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4. Mixing during halo mergers As structure in the Universe builds up through hierarchical 
clustering, dark matter halos merge and the metals contained in their hot gas compo- 
nents are mixed together. In this way, metals ejected by star-forming galaxies in the 
"field" at high redshift may end up contributing to the metals in the gas in rich clusters 
at 2; = 0. 

5. Mixing during galaxy mergers When two galaxies merge within a common dark matter 
halo, the metals contained in the cold gas and stars of both objects will be mixed in the 
merger remnant. We divide galaxy-galaxy merging events into two different classes: a) 
accretion events - this is when one galaxy is less than a third of the mass of the other 
one. b) major mergers - the galaxies are within a factor of three in mass. Major mergers 
are assumed to lead to the formation of elliptical galaxies. All the cold gas contained in 
the galaxies is turned into stars instantaneously in a "burst" and the enrichment during 
this phase is assumed to be a closed-box process, i.e. no metals escape from the galaxy. 
Accretion events do not induce any extra star formation. The metals present in the 
cold gas and stellar component of the satellite are simply added to the metals in the 
cold gas and stars of the central galaxy. Note that the factor of three is chosen so that 
the central galaxies of halos of circular velocity ~ 200 km s~^ will, on average, have 
bulge-to-disk ratios consistent with an Sbc-type spiral galaxy (see KWG). Our rules 
about star formation and enrichment during mergers and accretion events have been 
chosen for the sake of simplicity. They do not affect the results presented in this paper 
because, as we have mentioned, most stars in ellipticals do not form in these phases. 

To summarize, the main parameters of our model that control the metallicities of the 
stellar, cold gas and hot gas components in galaxies are the following: 

• Qft the density of the Universe in baryons 

• a the star formation efficiency 

• e the feedback efficiency 

• y the yield of heavy elements (instantaneous recycling assumed) 

• / the fraction of heavy elements ejected directly into the hot halo gas without first being 
mixed with the cold gas 

As in previous work (see KWG), we have adopted a "Milky Way" normalization, i.e. we 
tune the parameters so that the properties of the central galaxy in a halo of Vc = 220 km s~^ 
match those of our own Galaxy. In practice, the tuning procedure operates as follows. We first 
choose a value of and e control the luminosity and cold gas mass of the galaxy and are 

adjusted to obtain a S-band absolute magnitude of ~ —20.5 and a cold gas mass of 6 x IO^Mq. 
The yield is chosen so that the galaxy has a y-light weighted mean metallicity of 0.7 solar, a 
value consistent with the mean observed metallicity of stars in the solar neighbourhood (see 
Tinsley 1980 and references therein). It should be noted that in our model, since metals are 
driven out of a galaxy by winds and diluted in the surrounding gas, the mean metallicity of the 
stars in a galaxy will always be lower than the yield y. The higher the feedback efficiency, the 
higher the value of y required to obtain a metallicity of 0.7 solar. The parameter / controls the 
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total mass of metals contained in the hot halo gas. Since there are no reliable measurements 
of the metallicitics of galactic halo gas (indeed. X-ray emission from hot gas around spiral 
galaxies has yet to be detected), we tune / to match the abundance of metals observed in 
the intracluster medium. A compilation of data by Arnaud et al (1992) shows that the Fe 
abundance in a sample of nearby clusters is ~ 0.3 solar, with substantial scatter. We have 
thus chosen / so that the metallicity of the hot gas in halos of circular velocity 1000 km s~^ 
matches this value. Note that the properties of a Milky Way-type disk galaxies were explored 
in detail in a paper by Kauffmann (1996a) using a chemical evolution model very similar to 
the one described above. In that paper, it was demonstrated that the classic G-dwarf problem 
was solved by rapid early enrichment in the low mass objects that later merged to form the 
Galaxy, and that the predicted age-metallicity distribution of the stars in the Galactic disk 
agreed rather well with a recent compilation of data on solar neighbourhood stars. 

In Table 1, we list parameter values for two representative models. Model A (or the strong 
feedback model) has Qh = 0.1. For our adopted Hubble Constant of 50 km s~^ Mpc~^, this 
value lies at the upper end of the permitted range of according to current constraints from 
models of Big Bang nucleosynthesis (see Turner 1996 for a review), but at the lower end of 
measurements the fraction of baryons in rich clusters (White et al 1993; White & Fabian 
1995). Model A also has a high supernova feedback efficiency. The value of the yield (twice 
solar) is rather high, but still roughly consistent with theoretical calculations of the mass of 
different types of heavy elements ejected during Type II and Type la supernova explosions 
assuming a standard Salpeter IMF (Tsujimoto et al 1995 ). Model B (or the weak feedback 
model) has fib — 0.06. It has a considerably smaller feedback efficiency and yield. 

3 The New Spectral Synthesis Models 

We compute the spectrophotometric properties of model galaxies using new population synthe- 
sis models by Bruzual & Chariot (1997, in preparation). These span the range of metallicitics 
5 X 10"^ < Z/Zq < 5 and include all phases of stellar evolution from the zero-age main 
sequence to supernova explosions for progenitors more massive than 8Mq, or the end of the 
white dwarf cooling sequence for less massive progenitors. The models are based on recent 
stellar evolutionary tracks computed by Alongi et al. (1993), Bressan et al. (1993), Fagotto et 
al. (1994a, b, c), and Girardi et al. (1996), supplemented with prescriptions for upper-AGB 
and post-AGB evolution. The radiative opacities are taken from Iglesias et al. (1992). In the 
version used here, we adopt the library of synthetic stellar spectra compiled by Lejeune et al. 
(1997b, in preparation; see also 1997a) for all metallicitics. This library is based on spectra 
by Kurucz (1995, private communication; see also Kurucz 1992) for the hotter (0-K) stars, 
Bessell et al. (1989, 1991) and Fluks et al. (1994) for M giants, and AUard & Hauschildt 
(1995) for M dwarfs. The Lejeune et al. spectral library also includes semi-empirical correc- 
tions for blanketing, a well-known limitation of synthetic spectra (see, for example. Chariot, 
Worthey & Bressan 1996, and references therein). The resulting model spectra computed 
for stellar populations of various ages and metallicities have been checked against observed 
spectra of star clusters and galaxies (Bruzual & Chariot 1997; Bruzual et al. 1997). 

In addition, the models predict the strengths of 21 stellar absorption features, including 
the Mg2 index near 5000 A that is often used as a spectral diagnostic in early-type galaxies. 
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The predictions are computed using the Worthey et al. (1994) analytic fitting functions for 
index strength as a function of stellar temperature, gravity and metallicity that are calibrated 
empirically on a sample of 460 Galactic stars. This constitutes the standard "Lick/IDS" 
system, in which each index is defined by a central bandpass bracketted by two pseudo- 
continuum bandpasses. The integrated index of a stellar population is then computed by 
weighting the contributions from individual stars by their level of continuum (see Worthey et 
al. 1994, and Bressan et al. 1996 for more detail). It is worth noting that the most massive 
elliptical galaxies exhibit [Mg/Fe] ratios in excess of that found in the most metal-rich stars in 
the solar neighborhood (by ~ 0.2 — 0.3 dex; see Worthey, Faber, & Gonzalez 1992). While this 
may limit the accuracy of the predicted Mg2 indices of bright elliptical galaxies, the recent 
models of Bressan et al. (1997, in preparation) show convincingly that an enhancement in 
light elements at fixed total metallicity has virtually no effect on the other spectrophotometric 
properties of model stellar populations. 

A complete discussion of the differences between the spectrophotometric predictions of 
these models with those in previous studies will be presented in Bruzual & Chariot (1997). 
The typical discrepancies between the properties of stellar populations of the same input 
age and metallicity that are obtained by using the spectral synthesis models constructed by 
different groups of scientists, have already been illustrated by Chariot, Worthey & Bressan 
(1996). These can reach up to 0.05 mag in rest-frame B — V, 0.25 mag in rest- frame V — K 
and a 25% dispersion in the T^-band mass-to-light ratio. With these uncertainties in mind, we 
will concentrate more on understanding the trends seen in the observations than on obtaining 
exact fits to the data. 

4 Theory of the origin of a mass-metaUicity relation for 
eUipticals 

As was discussed previously, the assumption that supernova feedback ejects cold gas and 
metals from galaxies and that the efficiency of this process scales as where Vc is the 

circular velocity of the galaxy, means that more massive and luminous disk galaxies will be 
more metal rich. Because of the l/V^ scaling, the metaUicity-luminosity relation will saturate 
at large Vc when the feedback efficiency becomes neghgible and all gas and metals are retained 
in the galaxy. For large Vc, galaxies are effectively closed boxes. Thus in our scenario, metals 
in the hot gas component, including the gas in rich clusters, are produced primarily by low- 
mass galaxies since these lose their gas very easily. In clusters, most of the enrichment takes 
place at high redshift before the cluster has assembled (see section 8) . 

In figure 1, we show the metaUicity-luminosity relation for disk galaxies at 2; = in our 
model. We have plotted the l^-light weighted mean metallicities of the stars in these disks. 
The result for model A is shown as a solid line and that for B as a dashed line. As expected 
the relation rises much more steeply for model A, which has much more efficient feedback. 
In model B, the relation saturates at ^-absolute magnitudes of around —20. The model A 
relation is in better qualitative agreement with the observed metaUicity-luminosity relation 
of spiral galaxies (measured using emission lines from the gas in HII regions) since this shows 
no turnover (Roberts & Haynes 1994). 

In the merger model, it is not a priori clear that ellipticals should also obey a metallicity- 
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Figure 1: The metallicity-luminosity relation for disk galaxies in the model: the V-luminosity 
weighted mean metallicity of the stars is plotted as a function of V- band absolute magnitude. 
The sohd hne is for the high feedback model A and the dashed hne is for the low feedback 
model B. 

luminosity relation. For this to be the case, bright ellipticals have to be formed from the 
mergers of systematically bigger disks than faint ellipticals. Figure 2 demonstrates that this 
is indeed the case. We plot the mean progenitor mass of ellipticals in present-day clusters 
of circular velocity 1000 km versus their stellar mass. The mean progenitor mass is 
computed by weighting each star in the galaxy by the mass of the disk-galaxy progenitor in 
which it formed and is expressed in units of the total mass of the elliptical. As can be seen, 
the progenitors of elliptical galaxies of 2 x IO^^Mq are typically 10 times more massive than 
those of IO^^Mq ellipticals. The error bars show the 5th to 95th percentile spread in mean 
progenitor mass among ellipticals in each mass bin. 

In figure 3 , we show the metallicity-luminosity and age-luminosity relations for present- 
day cluster ellipticals. The left panel is for model A and the right panel is for model B. We 
plot P^-light weighted stellar quantities as before. Note that the ages differ somewhat from 
those given in Kauffmann (1996b), because the V-light weighting has changed as a result of 
the implementation of the new spectral synthesis models. The rms scatter in age is 1.2 Gyear 
and the scatter in metallicity at fixed luminosity is 0.14 in solar units. Model A again has 
a much steeper metallicity-luminosity relation than model B. Both models show very little 
trend in age with luminosity. If anything, the brighter ellipticals have somewhat younger ages. 
We conclude, therefore, that in the hierarchical merger model, the observed colour-magnitude 
relation for elliptical galaxies must arise as a metallicity effect. 
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Figure 2: The mean progenitor mass of ellipticals in clusters of circular velocity 1000 km s^^ 
is plotted as a function of their stellar mass. The progenitor mass is expressed as a fraction 
of the final mass of the eUiptical. Errorbars show the 5th-95th percentile range in progenitor 
mass. 

5 The colour-magnitude and Mg2-a relation for cluster 
ellipticals 

In figure 4, we show the U — V and V — K colour-magnitude relations for cluster ellipticals. 
Results for model A are shown in the left-hand panel and results for model B in the right. 
The straight lines are the fits to the U — V and V — K colour-magnitude relations of eUiptical 
galaxies in the Virgo and Coma clusters derived by Bower, Lucey & Ellis (1992, hereafter 
BLE). Model A fits the observations very well. We have computed least-squares estimates of 
the slopes and scatter of the Model A relations and compare them in Table 2 with the values 
derived by BLE for the Coma and Virgo ellipticals. Note that the galaxies plotted in figure 
4 are from 5 different Monte Carlo realizations of the formation of a cluster of = 1000 km 
s~^. Even so, the relation has very small scatter. Since model A fits the data so much better 
than model B, in future we show results only for model A. 

It is well known that the Mg2 line index correlates strongly with the stellar velocity disper- 
sion of elliptical galaxies (e.g. Bender, Burstein & Faber 1993). In our models, this relation 
is less straightforward to derive as we do not have a dynamical model for computing a. For 
simplicity, we have assumed a constant ratio Vc/a, where Vc is the circular velocity of the 
halo in which the elliptical formed during its last major merger. If the stars were a non-self- 
gravitating test population with density profile oc within a singular isothermal halo, one 
would have Vc/cr = 3^/^ (Gunn 1982). 

The model Mg2 indices, computed using the population synthesis models as described in 
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Figure 3: The metallicity-luminosity and age-luminosity relations of elliptical galaxies in 
clusters of circular velocity 1000 km at the present day . l^-luminosity weighted quantities 
are plotted, as before. 
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Figure 4: The U — V and V — K colour- magnitude relations of ellipticals galaxies in clusters 
of circular velocity 1000 km s~^. The sohd lines are the fits to the data on Coma and Virgo 
ellipticals taken from Bower, Lucey & Ellis (1992). 
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section 3, are global indices, i.e. they are averaged over the entire stellar mass of the galaxy. 
For nearby ellipticals, however, the observed indices come from measurements in the central 
regions of the galaxies. Since the abundance of magnesium is known to fall off quite strongly 
with radius, we are forced to correct for this effect if we want to make a meaningful comparison 
with the data. We have used the observed Mg2-radius relation derived by Davies, Sadler & 
Peletier (1993) (Mg2 oc — 0.059 log(it!/it!e), Re is the effective radius of the eUiptical), to 
transform the global indices to central values. Since our model does not make any prediction 
for the sizes of elliptical galaxies, we use the observed relation between ^-luminosity and 
effective radius given in Guzman, Lucey & Bower (1993) to assign a value of i?e to each of 
the model galaxies. 

The uncorrected and corrected Mg2 - a relations are shown in figure 5. The solid line is 
a fit to the data of J0rgensen, Pranx & Kjaergaard (1996). These authors quote Mg2 indices 
corrected to a circular aperture with diameter 1.19 h^^ kpc, equivalent to 3.4 arcsec at the 
distance of the Coma cluster. Since there is substantial scatter both in the observed Mg2 
abundance gradients of ellipticals and in the values of Re at fixed luminosity, which we ignore 
in making the corrections, we will not attempt to analyze the scatter in the relation. What 
is clear, however, is that the slope of the corrected relation matches the observations rather 
well. To fit the observations, we require Vc/cr ~ 2. 

In figure 6, we plot the Mb — cr Faber- Jackson relation for ellipticals (also assuming VJ a — 
2) . The line is a fit to the Mb — cr relation of Coma and Virgo cluster eUipticals from Dressier 
et al (1987). It is reassuring that we are able to fit both the Mb — cr and the Mg2 — a relation 
with the same parameters. This gives us confidence that the different outputs of the spectral 
synthesis models are giving consistent results. 

6 The mass-to-light ratios of cluster ellipticals 

If the luminosity profiles and the dynamical structure of all elliptical galaxies are similar, the 
virial theorem imphes that the mass-to-light ratio M/ L is a simple function of the effective 
radius Rg, the mean surface brightness (/)e, and the central velocity dispersion a. All three 
quantities may be easily determined observationally. 

It should be noted that the M/L ratio is a measure of both the stellar population and the 
dark matter content of galaxies. Our models can be used to predict the stellar M/L ratios of 
ellipticals. In figure 7, we show {M/L)v and [M/ L)k as a function of stellar mass. The lines 
are fits to recent data on Coma cluster eUipticals taken from Mobasher et al (1997). It should 
be noted that the values we obtain for our stellar M/ L ratios are very sensitive to the choice 
of IMF. The results shown in figure 6 are for a Scalo (1986) IMF with a lower- mass cutoff of 
0.1 Mq. a Salpeter IMF with the same lower cutoff gives M/L ratios a factor of 2-3 higher. 
The slope of the relation is, however, independent of the IMF. As can be seen, our models fail 
to reproduce the observed increase in M/L with mass. Our {M/L)k ratios actually decrease 
with mass, whereas Mobasher et al find {M/L)x oc M", with a — 0.14 ± 0.01. 

One way of resolving this problem would be for more massive ellipticals to contain a larger 
fraction of dark matter than less massive ellipticals. In our model, this requires that more 
effective mixing of dark matter with stars take place during the formation of more massive 
ellipticals. As can be inferred from figure 2, massive ellipticals have typically more merging 
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Figure 5: The corrected and uncorrected Mg2 — a relations computed from our models. In 
the lower panel, the Mg2 indices have been transformed from global indices to central indices 
as described in the text. The line is a fit to the data from J0rgensen, Pranx & Kjaergaard 
(1995). 
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Figure 6: The S-band "Faber- Jackson" relation for ellipticals in the model. The line is a fit 
to the data for Coma Cluster eUipticals taken from Dressier et al (1987). 
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Figure 7: The F-band and X-band mass-to-light ratios of elliptical galaxies in clusters are 
plotted as a function of their stellar mass. The hues are fits to data from Mobasher et al 
(1997). 
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events than less massive objects, so this appears plausible. Observationally, evidence for dark 
halos around ellipticals has been found through the studies of HI kinematics (e.g. Franx et 
al 1994), X-ray emission (e.g. Forman et al 1994), radial velocities of planetary nebulae and 
globular clusters (e.g. Mould et al 1990, Hui et al 1994), gravitational lensing (e.g. Maoz & 
Rix 1993) and measurements of the shape of the stellar hne-of sight velocity distribution (e.g. 
CaroUo et al 1995). More recently, Rix et al (1997) have analyzed the velocity profile of the 
elliptical galaxy NGC 2434 and show that roughly half the mass within an effective radius is 
dark. 

It should also be noted that recent analyses of the light profiles of elliptical galaxies show 
a shght breaking of the generally assumed homology. Allowing for this effect can substantially 
reduce the tilt of the M/L relation (Hjorth & Madsen 1996; Graham & CoUess 1997). This 
would also bring our results more in line with the data, and the variation in merger number 
with luminosity might offer an explanation for the variation in profile shape. 

Finally, Chiosi et al (1997) claim to obtain a good fit to the observed M/L relation using 
an IMF that varies as a function of the temperature, density and velocity dispersion of the 
gas that forms the eUiptical. Their model does not attempt to place eUiptical formation in its 
cosmological context. 

7 Evolution at high redshift 

In figure 8, we show age-luminosity and age-metallicity relations for elliptical galaxies in 
clusters of circular velocity 1000 km s"^ aX z = 0.4 and z = 1.5. As can be seen, the metallicity- 
luminosity relation does not evolve very much with redshift. By 2; = 1.5, it has somewhat 
larger scatter, but an My — —21 elliptical still has a metallicity of 0.6 solar on average, 
compared with 0.75 solar at the present day. The age-luminosity relation shifts downwards 
with redshift, and the magnitude of the shift is roughly to — t{z), where t{z) is the age of 
the Universe at redshift z. Note that the scatter in age is smaller at high redshift, reflecting 
the compression of the time interval over which the stars in the elliptical are constrained to 
form. Thus in our models, eUiptical galaxies of a given mass in high redshift clusters form 
and are enriched in much the same way as elliptical galaxies of the same mass in clusters at 
the present day. By selecting ellipticals in the richest chisters, one is automatically biasing 
ones sample of galaxies to the objects that turned around, collapsed and merged at the very 
highest redshifts. It is for this reason that the properties of these galaxies appear to follow 
the so-called "passive evolution" prediction (e.g. Ar agon- Salamanca et al 1993; Stanford et 
al 1995; Dickinson 1995; Franx & van Dokkum 1996; Bender, Ziegler & Bruzual 1996; Ellis 
et al 1996). In practice, however, the ellipticals in the high z clusters are not necessarily the 
direct progenitors of the Es in present-day clusters, so it is incorrect to set constraints on the 
formation of the present-day elliptical galaxy population using these objects (Chariot Sz Silk 
1994). This can only be done using a complete redshift survey of these objects (Kauffmann, 
Chariot & White 1996). 

In figure 9 we show the predicted restframe U — V colour-magnitude relation at four 
different redshifts. For reference, the solid line is the observed relation of BLE at z = 0. 
The dashed line indicates the shift in colour expected for a passively evolving population of 
solar metaUicity with a formation redshift of 10. Note that there is a progressive flattening 
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Figure 8: The metallicity-luminosity and age-luminosity relations of elliptical galaxies in 
clusters of circular velocity 1000 km s~^ at 2; = 0.4, z — 0.7 and z — 1.5 
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in the slope of the colour-magnitude relation with redshift, although it only becomes really 
noticeable at 2; > 1. This is caused by the fact that the bigger ellipticals have slightly younger 
mean ages in our model, and the effect of this on observables such as colour and magnitude, 
is much more pronounced at early times. The Mg2 — a relation also shows a flattening at 
z > 1, although it is a smaller effect than in the U — V colours. The slope and scatter of 
the rest-frame U — V colour- magnitude relation at each redshift are hsted in Table 2. The 
scatter in the rest-frame U — V colours of ellipticals in clusters at 2; ~ 0.54 has been analyzed 
in detail by Ellis et al (1997). They find no change with respect to the z — relation, which 
is in good agreement with our results. 

8 Metals in the intracluster gas 

Since our model follows the exchange of metals between the stars, cold gas and hot gas 
component of galaxies, we are in a position to study the origin of metals in the intracluster 
medium. Recall that we chose the parameter / to obtain a mean metallicity of 0.3 solar 
in the hot gas in a halo with circular velocity 1000 km s^^. Rather than the mass fraction 
of metals, many analyses have focused on the mass-to-light ratio of heavy elements in the 
ICM, as this quantity is independent of the degree to which the metals have been mixed 
with primordial gas. Most studies concentrate on iron, as this is the only element for which 
it has been possible, until recently, to derive accurate abundances using X-ray spectroscopy. 
For example, Renzini et al (1995) define the iron mass-to-light ratio (IMLR) as the ratio 
of the iron mass in the ICM to the total 5-band light of the cluster. They derive values 
of Mjf^^ /Lb in the range 0.01-0.02 Mq/Lq. The clusters we have modelled, which have 
mass 5 x lO^^M© and a baryon fraction of 10%, typically contain a total S-band luminosity 
of 2 X of which 60-70% is in early-type galaxies. Our total 5-luminosity and E/SO 

fraction agree well with the values measured for the Virgo cluster (Kraan-Korteweg 1981). We 
obtain M^^^ / Lb = 0.015 Mq/Lq, in good agreement with the observations. Recall that our 
model has a yield of twice solar. This is larger than the conventional value adopted in most 
analyses (eg Arnaud et al 1992). However, as discussed by Renzini et al, the uncertainties in 
theoretical estimates of the yield of nucleosynthesis products from supernova explosions are 
rather large. Recent calculations by Tsujimoto et al (1995) actually give a (Sn Il+Sn la) Fe 
yield close to twice solar, assuming a Salpeter IMF and the Sn la rate that best reproduces 
the abundance pattern of heavy elements in the solar neighborhood. We note that, according 
to Loewenstein & Mushotsky (1996), the abundances of elements such as O, Ne, Si and S in 
nearby clusters may be even higher than that of Fe, although the errors in the measurements 
are still quite large. The theoretical interpretation of these elements is in principle much 
simpler than that of Fe, because they are produced almost exclusively in Sn II explosions. 
Matching these results might require yields even higher than twice solar and would most 
hkely only be compatible with a top-heavy IMF. We leave a more detailed examination of 
these issues to future work. 

We can also use our models to study the history of the enrichment of the ICM and the 
distribution in mass (or circular velocity) of the galaxies that ejected the metals. This is 
illustrated in figure 10, where we plot the fraction of the total metals present in the present 
ICM that was ejected at a given redshift and from a galaxy of given circular velocity. As 
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Figure 9: The rest-frame [/ — colour-magnitude relations of elliptical galaxies in clusters 
of circular velocity 1000 km s~^ at z = 0.4, 0.7, 1.5 and 2.0. The solid lines are the z = 
data from Bower, Lucey & Ellis (1992). The dashed lines show the shift in colour expected 
for a passively evolving population with Zform = 10- 
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Figure 10: Top panel: the fraction of metals in the intracluster medium of a present-day 
cluster with circular velocity 1000 km s~^ ejected by disk galaxies with circular velocity T^- 
Bottom panel: the fraction of metals in the ICM today that were ejected at reshift z. 

can be seen, more than 80% of the metals were ejected at redshifts greater than 1, 60 % at 
redshifts greater than 2 and 40% at redshifts greater than 3. More than 40 % of the ejection 
was done by galaxies with circular velocities less than 125 km s~^, i.e. "dwarf" galaxies have 
contributed significantly to the metals in the ICM. Because the enrichment predominantly 
occurred at high redshift, our models predict that abundances in the ICM should evolve very 
little with redshift. Indeed, hy z = 1.5, the model predicts that the ICM metallicity is still 
~ 0.25 solar. Our results are in good agreement with the lack of evolution seen in the Fe 
abundances of clusters out to z ^ 0.3 as measured by the ASCA satellite (Mushotsky & 
Loewenstein 1997), although the model abundances show significantly less scatter than the 
data, suggesting that star formation in real cluster galaxies may be more stochastic than we 
have assumed. 

9 Discussion 

The failure to match the slope of the colour-magnitude relation of elliptical galaxies has 
been a serious defect of all previous attempts to model the formation of galaxies within the 
framework of a hierarchical cosmology (Lacey et al 1992 ; Kauffmann 1996b; Baugh, Cole & 
Frenk 1996). We have demonstrated that the inclusion of chemical evolution in the model 
solves this problem. With a physically realistic choice of parameters we now obtain good fits 
to the U — V and V — K colour-magnitude relations as well as the Mg2 — a relation. It is 



T ' ' ' r 



1 , 

J 1 1 1 L 




19 



interesting that the colour-magnitude relation only extends to the brightest ellipticals if we 
require that feedback be very efficient, even in relatively massive galaxies. Using the value of 
e given in Table 1 for model A, we can rewrite the cold gas reheating rate in equation 2 as 




280A;ms-^ 



(3) 



dt 



This implies that the reheating rate is 60% larger than the star formation rate for our own 
Milky Way. Because of the efficient feedback, our best fit model also requires a high value of 
Jib, so that enough gas is able to cool and form the stars of the brightest galaxies. 

Note that we have studied only global colours and metallicity indices values in this paper. 
It is well-known that elliptical galaxies exhibit strong gradients in their observed properties. 
In particular the nuclei of ellipticals host stellar populations with high metallicity ( ~ 0.5 to 
~ 3 Zq), a magnesium overabundance with respect to iron, and a small spread in metallicity. 
As discussed by Greggio (1997), the first two characteristics of the nuclei can be explained 
if gas is processed through multiple generations of stars on a short timescale. The small 
metallicity spread arises if the gas turning into stars in the nuclei has been pre-enriched. In 
N-body plus hydrodynamics simulations of the merger of two disk galaxies, gas is seen to 
be driven very effectively to the centre of the merger remnant, where it will likely turn into 
stars on a short timescale (Barnes & Hernquist 1996). Whether this picture can account for 
the observed properties of elliptical nuclei has yet to be investigated in detail. In principle, 
it would be an interesting test of our models, which are able to predict when the mergers 
occurred, and the metallicity and mass of the gas in the disk-galaxy progenitors. 

As in many other analyses, we fail to account for the tilt of the M/L relation, if M is 
purely a measure of the stellar mass of the galaxies (Djorgovski & Santiago 1993; Worthey 
1994). As we have discussed, there are several possible solutions to this problem, but one 
interesting one would be if more massive ellipticals contained a greater percentage of dark 
matter in their inner regions than less massive ellipticals. This might occur if dark matter 
is more efficiently mixed with the stars during the formation of the more massive systems. 
Clearly this is a problem that needs to be addressed using N-body simulations that take into 
account the full dynamics of the problem. 

We have also predicted that the mass-metallicity relation for elliptical galaxies in clusters 
remains virtually unchanged out to high redshift. It is only the age of the galaxies that ought 
to decrease. This is a prediction that might be testable using stellar population indicators that 
are able to disentangle age and metaUicity effects (eg Jones & Worthey 1995). In addition, 
we predict that there should be a slight fiattening in the slope of the colour-magnitude and 
Mg2 — o" relations at redshifts greater than 1, reflecting the fact that more massive ellipticals 
are slightly younger in age than the less massive ones. 

Finally, X-ray satellites, such as ASCA, are beginning to provide a wealth of information 
about the abundances many different elements in the intracluster gas, as well as about the 
evolution in ICM abundances to high redshift (e.g. Mushotsky & Loewenstein 1997). This 
information will provide further constraints on the models, in particular on our prescriptions 
for star formation and supernova feedback. These are indeed the most uncertain elements of 
our galaxy formation scheme. 
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Table 1: The parameters for the models 



a e y f 
Model A 0.1 0.2 0.4 2.0 0.3 
Model B 0.06 0.1 0.05 1.2 0.45 



Table 2: Best-fit regression parameters 
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